clear; clc; close all;

% 参数设置
q = 0.4;
a = 2.3617;
b = 1.2;
P = 1;
Q = 1;
N = 5000; % 迭代步数，足够长展示吸引子形态

% 两组初始条件
x0_1 = -0.5; y0_1 = 0.5;
x0_2 = 3;    y0_2 = -3;

% 计算轨迹
[x1, y1] = SCLMM(q, a, b, P, Q, x0_1, y0_1, N);
[x2, y2] = SCLMM(q, a, b, P, Q, x0_2, y0_2, N);

% 绘图
figure; hold on;
plot(x1, y1, '.', 'MarkerSize', 1, 'Color', [0 0 1]); % 蓝色点，初始值1轨迹
plot(x2, y2, '.', 'MarkerSize', 1, 'Color', [1 0 0]); % 红色点，初始值2轨迹

xlabel('x(n)', 'FontSize', 14);
ylabel('y(n)', 'FontSize', 14);
title('Fig.8(a) 2D-SCLMM 多稳态共存吸引子', 'FontSize', 16);
grid on; axis equal;
box on;
